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The CDMS and EDELWEISS collaborations have combined the results of their direct searches 
for dark matter using cryogenic germanium detectors. The total data set represents 614 kg -days 
equivalent exposure. A straightforward method of combination was chosen for its simplicity before 
data were exchanged between experiments. The results are interpreted in terms of limits on spin- 
independent WIMP-nucleon cross-section. For a WIMP mass of 90 GeV/c 2 , where this analysis is 
most sensitive, a cross-section of 3.3 x 10 -44 cm 2 is excluded at 90% CL. At higher WIMP masses, 
the combination improves the individual limits, by a factor 1.6 above 700 GeV/c 2 . Alternative 
methods of combining the data provide stronger constraints for some ranges of WIMP masses and 
weaker constraints for others. 
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The problem of dark matter has been an open question 
since 1933 [lj : most of the matter in the Universe appears 
only through its gravitational interactions. Evidence sug- 
gests that this dark matter may be made up of weakly 
interacting, massive particles (WIMPs) [2]. Supersym- 
metric theories, and other extensions of the standard 
model of particle physics, predict plausible candidates 
for WIMPs 0, 0] , and experimental efforts have been un- 
der way since the mid-1980s to detect them 0,0, @. The 
challenge is great because the average deposited energies 
involved are quite low, between a few keV and a few tens 
of keV, and because the event rates are minute compared 
to normal levels of radioactive backgrounds. Therefore, 
experiments searching for WIMPs are typically located 
underground to reduce exposure to cosmic radiation, and 
they generally deploy detectors with some form of back- 
ground identification. For instance, cryogenic detectors 
operating at tens of millikelvins use a simultaneous mea- 
surement of phonons and charge to efficiently reject most 
of the dominant radioactive background (e.g. gamma 
particles) which interacts with electrons and, for a given 
energy deposit, ionizes more than the nuclear recoils in- 
duced by WIMP scattering. The experiments CDMS 0- 
G3 and EDELWEISS have set strong constraints 

on WIMPs over the past decade using this technique with 
germanium targets. Other experiments use other target 
nuclei and techniques fl4T-[l7T] . Use of the same target nu- 
cleus and experimental technique by CDMS and EDEL- 
WEISS offers the possibility of combining their results 
to establish stronger constraints on WIMPs without in- 
troducing more model dependence than the individual 
results already had. 

Published results @-G3 from the CDMS experiment 
were obtained at the Soudan Underground Laboratory 
in Minnesota, a site of intermediate depth at 2100 me- 
ters water equivalent (mwe). CDMS has operated be- 
tween six and thirty cryogenic semi-conductor detectors. 
In addition to the dual phonon-ionization measurement 
used to reject the majority of the highly-ionizing electron 
recoil background in the bulk of the detectors, CDMS de- 
vices measure non-thermal phonon information to iden- 
tify background near the surface of the detector that 
could otherwise be confused with signal [l8|]. Moreover, 
all installed detectors have been used to identify back- 
grounds through multiple scattering. However, only the 
germanium ones (each of mass ~ 230 g) that were fully 
functioning (between four and fifteen at various times) 
were used to set constraints on the WIMP cross section. 
Thresholds as low as 5 keV have been obtained in im- 
proved analyses [l9[ of certain runs 0, [l(| ■ The effective 
net exposure of the CDMS data set as a function of recoil 
energy, obtained by multiplying the energy-de pen dent 
efficiency by the exposure, is shown in Fig. [T] [2(|. It 
reaches a maximum of 379 kg-d at 25 keV. It then drops 
by roughly a quarter to the value obtained at 100 keV, the 



highest energy considered in these CDMS searches. To 
avoid bias, these cuts were set in a blind manner, without 
knowledge of events in, or near, the signal region. The 
four remaining candidate events observed in the CDMS 
data are all below an energy of 20 keV (see Table U for 
energies of individual events). The expected number of 
background events misidentified as signal candidates in 
the full CDMS data set is ~ 2 events [jj, 1§J, [l9| , with sig- 
nificant systematic uncertainty on the background at the 
lowest energies [19j |. 
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FIG. 1. Top: effective net exposure curves for CDMS (dotted 
red curve), EDELWEISS (dash-dotted blue curve), and for 
the combined experiments (dashed black curve). Middle: en- 
ergies of the five observed candidate events in EDELWEISS 
(blue lines). Bottom: energies of the four observed candidate 
events in CDMS (red lines). All are below the EDELWEISS 
threshold. In the simple method of merging the two sets of 
data, experimental provenance of the individual events is not 
kept. 



TABLE I. Energies of the candidate events observed in each 
experiment [2C( ]. 



CDMS 



Event energy [keV] 



EDELWEISS 



5.3 7.3 12.3 15.5 



20.8 21.1 21.8 23.2 172 



The EDELWEISS experiment [11| has deployed ten 
germanium detectors of 400 g at the Modane Under- 
ground Laboratory, a deeper site with a rock overburden 
of 4800 mwe. These detectors also use phonons and ion- 
ization to reject bulk background. In contrast to CDMS, 
the EDELWEISS detectors use patterned electrodes cre- 
ating an inhomogeneous electrical field to identify surface 



background 24|. The corresponding effective net expq- 
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sure as a function of recoil energy is shown in Fig. [T] 20] . 
The distribution is very close to a step function with a 
sharp threshold of 20 keV and a flat plateau at 384 kg-d 
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90% CL Limits: Simple Merger of CDMS and EDELWEISS Data 
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FIG. 2. Top: 90% CL optimum interval upper limits on 
spin-independent WIMP couplings to nucleons as a function 
of WIMP mass, from the individual CDMS (red dashes) 
and EDELWEISS [IU (EDW, blue crosses) experiments, and 
from their simple merger (continuous black line). Also rep- 
resented are limits from the XENON 100 [lj (brown boxes), 
XENON 10 [U (green crosses V.CRESST II pj] (brown dot- 
dashed line) and ZEPLIN III [jLy] (pink dots) experimentjq, 
and supersymmetric parameter-space predictions [22l [23 | 
(filled gray regions). Bottom: gain obtained from the straight- 
forward combination with respect to the strongest individual 
limit of CDMS and EDELWEISS (in effect that of CDMS). 
Below masses of 50 GeV/c 2 , the combined limit is weaker 
than the best individual one; at higher masses, the gain is up 
to a factor 1.57. 

a The limits for XENON 10, CRESST II and ZEPLIN III use 
slightly larger halo escape velocities than this work 
(600-650 km/s). 



up to 200 keV. Though not blind, cuts for the EDEL- 
WEISS data were set without considering events in the 
signal region. Four of the candidate events observed in 
the EDELWEISS data set have energies in the 20-25 keV 
range; the fifth event lies at 172 keV. Precise values of 
energies are listed in Table |U The expected contribu- 
tion of known backgrounds to these events is at most 
3 events [ll|. Since the statistical sample is small and 
uncertainties in the backgrounds are still under investi- 
gation, all events are considered WIMP candidates for 
the purpose of deriving an upper limit on WIMPs. Lim- 
its obtained by EDELWEISS are very similar to those 
of CDMS at high WIMP masses, though weaker at low 
WIMP masses. 

Despite similar maximal exposures, the zero- 
background sensitivities of CDMS and EDELWEISS, 
estimated from the exposure- weighted efficiencies alone, 
differ significantly due to the higher threshold of the 
latter. For example, the 10-100 keV spectrum-averaged 
equivalent exposure for a WIMP of mass 90 GeV/c 2 
is 361 kg-d for CDMS and 253 kg-d for EDELWEISS. 



The corresponding values calculated for a mass of 
1 TeV/c 2 are 357 kg-d and 309 kg-d, respectively. 
They converge as the recoil energy spectrum hardens, 
explaining in large part the similar limits obtained 
at high mass. The rest of the difference between the 
limits obtained by the two experiments is due to the 
number and energies of the nuclear-recoil candidates 
observed in each. Previously published limits from each 
experiment were calculated individually using Yellin's 
optimum interval method [25[, which effectively sets a 
limit using the region of the observed spectrum that is 
most constraining for the given model, and applies the 
appropriate statistical penalty for choosing this region. 
All events are considered as potentially valid WIMP 
candidates; no expected background is subtracted. This 
technique has the advantage that it does not require 
any a priori information on the backgrounds, and is well 
suited to cases with significant systematic uncertainties. 

There are many possible procedures for combining re- 
sults from experiments based on the optimum interval 
method [26j . Without some assumption about the back- 
grounds in the two experiments, it is impossible to de- 
termine which method has the best expected sensitivity. 
Because not all background estimates for the full datasets 
considered had been made "blindly" before the signal re- 
gions were examined, and because of the uncertainties 
associated with these estimates, the collaborations de- 
cided not to try determining which method has the best 
expected sensitivity, but instead to use the same method 
each individual experiment had previously used 0, [H| 
for combining data from different detectors and differ- 
ent data sets. In this simple merging method (26|, the 
exposure-weighted efficiencies of both experiments are 
simply summed, and the events all treated on equal foot- 
ing without consideration of their experiment of origin. 
The experimental limits are then derived from the com- 
bined data set using the optimum interval method. The 
decision to use this method was made before data were 
exchanged, and indeed before the final analysis of the 
EDELWEISS data was available. 

The combined effective net exposure curve is shown 
in Fig. □ The derived 90% CL limit on the spin- 
independent cross-section for WIMP-nucleon scattering, 
as a function of mass, is shown in Fig. [2] We have 
assumed a standard astrophysical halo model (WIMP 
mass density 0.3 GeV/cm 3 , most probable WIMP veloc- 
ity with respect to the galaxy 220 km/s, mean velocity 
of Earth with respect to the galaxy 232 km/s, galactic 
escape velocity 544 km/s [13]), and the Helm form fac- 
tor . The limit is strongest at a mass of 90 GeV/c 2 
with a value of 3.3 x 10~ 44 cm 2 ; it increases by less than 
4% over the 75 to 115 GeV/c 2 range. Figure [2] also com- 
pares this limit with the results from the individual ex- 
periments. The combined limits are stronger than those 
from the more sensitive of the two individual experiments 
(CDMS) for WIMP masses greater than 50 GeV/c 2 . The 
gain is as much as a factor 1.57 at the highest masses, 
compared to a maximum allowed gain of 1.9 from the 
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relative experimental exposures alone. In the optimum 
interval technique, the gain results from the large event- 
less interval between 23.2 and 172 keV. Below a mass 
of 35 GeV/c , the average recoil energy expected for a 
WIMP signal falls below 12 keV, and events start to ap- 
pear in the optimum interval, further degrading the limit. 
At high masses, the optimum interval contains no events, 
and the abrupt change in the limit is due to subtleties of 
the optimum interval method. Overall, the upper limits 
correspond to an expected number of detected WIMPs 
ranging from roughly 12 at low masses to about 4 at high 
masses. 

While this combination procedure does not preserve 
the identification of the candidate by its experiment of 
origin, and thus cannot take account of differing back- 
grounds, other methods can. Such methods could be of 
interest in future combinations of data sets with signifi- 
cantly different expected backgrounds, if those are known 
a priori. For comparison, Fig. [3] shows the results of 
two alternate procedures [26| that may be expected to 
yield more constraining combined limits than the simple 
method used here. 

The first alternate procedure (the "minimum limit" 
method [26| ) in essence computes the upper limit from 
each experiment separately, takes the lower of the two 
limits, and applies the appropriate statistical penalty for 
the choice. This method is suitable if both experiments 
are background-limited. Figure [3] illustrates that this 
method results in less sensitive limits for the case at hand 
for high WIMP masses, since it effectively uses the results 
of only one experiment even for masses for which both ef- 
fectively have negligible candidate events. For low WIMP 
masses, however, it provides more stringent limits than 
the simple method. This result also is to be expected: 
for masses for which the number of candidate events is 
large for both experiments, this method effectively deter- 
mines which experiment provides the better limits and in 
practice ignores the other. 

The second alternate procedure (the "probability prod- 
uct" method [26J) considers, at each WIMP mass, the 
product of the probabilities that each experiment ex- 
cludes a given cross section as too high. This method has 
poor sensitivity if experiments are background-limited, 
since an experiment with higher background will be 
weighted equally, worsening the limit of the more sensi- 
tive experiment. It may have improved sensitivity com- 
pared to the simple method if both experiments are 
exposure-limited, yet have differing backgrounds. As il- 
lustrated in Fig. |31 this method gives a weaker limit than 
simple merging at low masses. However, as Fig. [3] also 
shows, for masses above 25 GeV/c 2 , this method pro- 
vides more stringent limits than the simple merger. This 
result occurs because the optimum intervals for CDMS 
and EDELWEISS arc different at these masses, so the ex- 
tra freedom (compared to the simple method) to choose 
different intervals for the two experiments results in a 
greater sensitivity. We note lastly that this procedure is 
not defined for WIMPs too light for at least one of the 



90% CL Limits: Alternative Combinations of CDMS and EDELWEISS Data 
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FIG. 3. Comparison of combined CDMS-EDELWEISS lim- 
its obtained by different statistical methods. The "mini- 
mum limit" method (orange crosses) is more constraining 
than the adopted, "simple merging" method (black curve) 
for low WIMP masses, but less constraining at high masses. 
The "probability product" method (purple circles) yields the 
strongest limits for masses above 30 GeV/c 2 , but has the 
weakest limits for low masses, and is not defined for WIMPs 
with masses below ~ 12 GeV/c 2 since such WIMPs cannot be 
observed in EDELWEISS due to its the higher threshold. The 
limit of the combined experiments under the simple method 
had no events been detected (red dashes) lies below these 
limits by a factor from ~ 1.5 to ~ 9. Bottom: gains of the 
alternative methods relative to the simple merger. 



experiments to observe, as is the case for EDELWEISS 
and WIMPs below ~ 12 GeV/c 2 . 

The difference between the result of the minimum limit 
procedure and that of the probability product calcula- 
tion suggests that the comparison of the two data sets 
could yield information on the backgrounds of each ex- 
periment, independent of that determined earlier (7L[Tl|. 
For each experiment, a likelihood was calculated as a 
function of WIMP mass using the experimental energy 
distribution of events and the expected distribution of 
WIMP events, and assuming no background. The EDEL- 
WEISS event at 172 keV is ignored since, even for a heavy 
10 TeV/c 2 WIMP, the probability that at least one of the 
five EDELWEISS events occurs at 172 keV or above is 
less than 2%. For CDMS and EDELWEISS individually, 
the maximum likelihood is obtained for WIMP masses 
below 17 GeV/c 2 ; however, for such masses, the expected 
rate in CDMS is much higher than that in EDELWEISS, 
in contradiction with the measurement. To further inves- 
tigate compatibility, we performed a likelihood ratio test. 
It rejects the hypothesis of no background at a confidence 
level greater than 99.8%. This result is fairly insensitive 
to changes in the WIMP astrophysical distribution, for 
instance increasing the escape velocity up to 650 km/s, 
or varying the average WIMP velocity between 150 and 
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350 km/s. Note however that no background subtraction 
has been performed to establish the limits in this work. 

In conclusion, we have improved constraints from sub- 
kelvin germanium detectors on WIMPs of mass greater 
than ~ 50 GeV/c 2 by a factor of up to ~ 1.6 thanks to 
a simple merger of the data from the CDMS and EDEL- 
WEISS experiments. Alternative methods that exploit 
the provenance of events could provide even stronger con- 
straints at certain masses. Except for simple merging, 
these methods apply without change to combining results 
from experiments with different target nuclei for a given 
set of astrophysical and physical assumptions. For simple 
merging of experiments with different recoil ener gy s cales 
due to different target nuclei, it is recommended [26| that 
before merging, the recoil energy of each event be trans- 



formed into its experiment's cumulative detection prob- 
ability function at that event's energy. 
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